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CI observations in the CQTau proto-planetary disk: 
evidence for a very low gas-to-dust ratio ? * 
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ABSTRACT 



Context. Gas and dust dissipation processes of proto-planetary disks are hardly known. Transition disks between Class II (proto- 
planetary disks) and Class III (debris disks) remain difficult to detect. 

Aims. We investigate the carbon chemistry of the peculiar CQ Tau gas disk. It is likely a transition disk because it exhibits weak CO 
emission with a relatively strong millimeter continuum, indicating that the disk might be currently dissipating its gas content. 
Methods. We used APEX to observe the two CI lines 3 />i -> 3 P at 492 GHz and 3 P 2 -> i P\ at 809 GHz in the disk orbiting CQ Tau. 
We compare the observations to several chemical model predictions. We focus our study on the influence of the stellar UV radiation 
shape and gas-to-dust ratio. 

Results. We did not detect the CI lines. However, our upper limits are deep enough to exclude high-CI models. The only available 
models compatible with our limits imply very low gas-to-dust ratio, of the order of a few, only. 

Conclusions. These observations strengthen the hypothesis that CQTau is likely a transition disk and suggest that gas disappears 
before dust. 
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1. Introduction 

Gas and dust disks orbiting Pre-Main-Sequence stars are be- 
lieved to be the birth-place of planetary systems. Several disks 
surrounding young stars have been detected and studied both 
in c ontinuum and mo lecular lines at mm/submm wavelengths 
(e.g. lPietu et al. 2007:, and references therein). However, little is 
known about how their gas and dust content dissipate. Transition 
disks between Class II (proto-planetary disks) and Class III (de- 
bris disks) remain hard to detect. Accretion, viscous spread- 
ing, planet formation, and photo-evaporation (enhanced by grain 
growth and dust settling) should play a role in the disappearance 
of the gas and dust disk, but their relative importance has not 
been estimated yet. 

Disks with strong mm continuum excess and weak CO lines 
are interesting in this respect, as they may represent a stage in 
which gas is being dissipated, while large dust grains respon- 
sible for the mm continuum still remain, perhaps settled in the 
mid-plane. The first case discovered was BP Tau which exhibits 
a warm (50 K), s mall (radius ~ 120 AU) disk with a very low 
CO to dust ratio dDutrev et alJl2003[). corresponding to a n ap- 
parent depletion factor 100 for CO. IChapillon et alJ d2008l here- 
after Paper I) have recently found that the Herbig Ae (HAe) stars 
CQTau and MWC758 have disks with similar characteristics. 
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IChapillon et all (l2008h showed that such a low CO to dust 
ratio can be explained by photo-dissociation by UV photons 
coming from the star with a standard gas-to-dust mass ratio. 
However, solutions with low gas-to-dust ratio also exist. These 
solutions differ by the amount of atomic and/or ionized carbon 
which result from the CO photo-dissociation. Thus, detection of 
either CI or C + provides a useful diagnostic on the disk struc- 
ture. An upper limit on the CI 3 P2 — > 3 Pi transition i n the Herbig 
Be sta r HD 100546 has been recently reported by iPanic et all 
(120101) . 

We report here a search for CI in CQ Tau. We present our ob- 
servations in Section 2. Our results are discussed and compared 
with chemical modeling in Section 3. We summarize in the last 
section. 



2. Observations and results 

We selected CQ Tau as our main target, because it appears to be 
the warmest source: the mm continuum and 12 CO observations 
of Paper I indicate a gas temperature > 50 K, a disk outer ra- 
dius near 250 AU, and an apparent CO depletion of 100. CQ Tau 
(see Table 1) is one of the oldest HAe star (~ 6 - lOMyr, al- 
though a revision of its distance implies that CQ Tau is younger 
than initially thought, see Paper I) surrounded by a resolved disk 
dMannings & Sargent|[T997t iDoucet et ail 120061). T his star ex- 
hibits an UX-Ori-like variability ( Natta et al. 1 19971) . in contra- 
diction with the observed low disk inclination (Paper I). Its dus t 
emissivity exponent is about 0.7 (Paper I and lTesti et al.ll2.Q03h , 
indicating significant grain growth. 
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Table 1. Properties of CQ Tau, also named HD 36910. 



Spec, type 


T,//(K) 


Stellar lum. ( L Q ) 


Distance'*' (pc) M/** 1 ( M Q ) 


Vta .'"» (km/s) 


PA ( ) ( " ) 


i cr*> R „ f (AU) ( **» 


A8-F2 


7200 


8-16 


140 1.8 


6.2 ± 0.04 


-36.7 ±0.3 


29.3 ±1.7 210 ±20 



w assumed to be the Taurus distance. ( ** ) derived according to interferometric 12 CO J=2-l data (Paper I) 
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> 3 Po (a) and 3 P2^> 3 P\ (b) spectra towards the 
Table 2. Sources coordinates (J2000). 



source 


RA 


DEC 


CQTau 


05:35:58.485 


24:44:54.19 


Mars w 


05:57:24 


23:30:41 


Orion-Bar 


05:35:15.09 


-05:26:36.6 



Aug.,25'" 2009 at OhOO (UT). 



We observed simultaneous ly the CI 3 Pi— > 3 Pq and 3 Pi^> 3 P\ 
lines in CQTau with FLASH (iHevminck et al.ll2006l) on APEX 
during August 2009 under excellent weather conditions (pwv 
< 0.3mm). We used wobbler mode (60"). The total integration 
time on source is 67min. Pointing and focus were frequently 
checked on Mars (the planet was very close to the source, see 
Table|2]i. Receiver tuning was checked on the Orion-Bar (Fig.[TJ. 
Our observation of Cl 3 Pi— > 3 Pq towards the Orion Bar peaks 
at T^ 17. 5 /if. Taking the beam efficiency measured toward 
the Moon ~ 0.8 (the emission being extended) this leads to a 
peak temperature of Tmb ~ 21 K. This is in quite good agree- 
ment with pre vious observations o f CI at 492 GHz toward the 
Orion Bar by iTauber et al] ( 1 1995b with the CSO 10.4m dish 
(15" beam, our pointing is between their "d" and "e" spectra, 
i.e. 17 < Tmb - 25 K, see their Figs. 1 and 2). 

The spectra towards CQTau are presented in Fig. [2] 
We did not detect any CI emission. The lcr sensitivity 
on the integrated area is given by the following formula 
cr(Jy.km.s _1 ) = S (Jy/K)<x(K) VtfvAv with <5v the spectral reso- 
lution and A v the line width. We take A v = 5 km.s -1 according 
to the PdBI CO observations (Paper I). The rms noise in Fig. [2] 
are (lcr) rms 0.020 and 0.046 K (Ta) for a spectral resolution of 
1.04 and 0.95 km.s 1 at 492 and 809 GHz, respectively. Taking 
the conversion factors between Jansky and Kelvin equal to 48 
and 70 Jy/K, we thus obtain 3<x upper limits of 6.6 Jy.km.s~ 1 and 
25 Jy.km.s~ 1 for the 3 P l ^ 3 P and 3 P 2 ^ 3 P\ lines (see Tableg). 



3. Discussion 

3.1. Chemical Modeling 




We u s e the PDR code from the Meudon group dLe Bourlot et al. 
U993HLePetit et alJl2006l) with the modification of the extinc- 



-40 -20 20 40 -40 -20 
Velocity (km/s) 

Fig. 2. (a) CI 3 Pi-> 3 P spectra, rms=0.020K and (b) CI 
3 P2-^> 3 P\ spectra, rms=0.046K towards CQTau. The spectral 
resolution is about 1 km/s. The dotted line correspond to the V/ Jr 
of the source. 

Table 3. Elemental abundances with respect to total hydrogen 



Element 


abundance 


C 


1.38 x 10"* 


O 


3.02 x 10~ 4 


N 


7.95 x 10- 5 


Mg 


1.00 x 10~ 8 


S 


2.00 x 10~ 6 


Si 


1.73 x 10~ 8 


Fe 


1.70 x 10- 9 



tion curve calculation according to the grain size distribution de- 
scribed in Paper I. The chemical network an d elemental abun- 
dances (given in Table [3) are similar to that of lGoicoechea et aD 
J2006I) . No freeze-out onto grains is considered, as the gas tem- 
perature determined by CO observation is > 50 K. We solve 
the radiative transfer perpendicular to the disk plane for each 
radii, with an incident UV flux scaling as 1/r 2 . Radiative trans- 
fer, chemistry and thermal balance are consistently calculated. 
For more details and a description of the use of the code see 
Paper I. 

From Paper I, our best models to reproduce the CO depletion 
and gas temperature in CQTau are 1) a disk illuminated by 
a strong UV field from the star (Draine field with a scaling 
factory = 10 4 at R= 100AU) with a normal gas-to-dust ratio 
(g/d = 100, model A) or 2) a disk illuminated by a weaker UV 
field {x = 10 2 ) with a modified g/d ratio of 10 (model B). In 
both cases we considered big grains with a maximum size of 
1 cm as suggested by the previous dust and CO observations. 
The computed surface densities and abundances of CO, CI and 
CII are shown in Figs. [3] and |4] left and middle panels. A cut 
perpendicular to the disk mid-plane in our models A and B 
shows three parts (FigHJi. At 100 AU radius, C + is the most 
abundance form of carbon in the upper layer, then there is a layer 
transition where CI is quasi-dominant and finally CO is the main 
carbon-component in the mid-plane. At larger radii (200 AU), 
C + stays however the dominant carbon-component even in the 
disk mid-plane, because the total opacity perpendicular to the 
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Fig. 3. Radial distribution of the surface density of C + , C, and CO for three models which are in agreement with the CO observations 



model A: g/d rotio=100 UV: x=10 4 model B: g/d rotio=10 UV: x=10 2 model C: g/d rotio=10 UV: Kurucz 
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10 1 1 10 1 1 10 1 1 

Av(mag) Av(mag) Av(mag) 

Fig. 4. Vertical distribution of the abundance of C + (blue), C (green), and CO (red) and gas temperature (pink dashed lines) for the 
three models which are in agreement with the CO observations. 



disk mid-plane is low. 



3.2. The UV problem 

As indicated in Paper I (§5.2), the actual UV spectra of HAe 
stars are not well known. Observations of MWC758 indicate a 
UV spectrum shortwards of 1500 A whose shape is better rep- 
resented by a Draine field than by a purely photospheric spec- 
trum (Paper I). From spectra obtained by the IUE satellite, the 
CQTau UV flux i s 10-30 times lower th an that of MWC758 
dGradv et alj|2005t iBlondel & Djiell2006l) , but the shape is un- 
known. To check the importance of the UV spectrum profile, 
we performed another model assuming a photospheric spectrum 
from a A3 type star (Kurucz database^). CQTau is a somewhat 
cooler star (A8 to F2), so this is most likely an overestimate of 
the UV flux. 

The CO abundance appears very sensitive to the UV profile. 
In absence of UV excess and in the case g/d=100 the chemistry 
is dramatically affected. The disk is mainly molecular and the 
C/CO transition occurs at low opacities (i.e. closer to the disk at- 
mosphere). This model leads however to too high CO surface 
densities (see Fig. [5] left panel). With g/d=10 (Model C, see 
Figs[3]-0, the CO surface densities are of the same order of 
magnitude as in model B and therefore also compatible with CO 
observations frome Paper I, but now the gas is warmer and there 
is more CI than C + (CI being the main carbon form in the mid- 
plane^ 

iBergin et al.l d2003l) have shown that the UV excess of 
T Tauri stars is dominated by strong lines emission, in particu- 
lar Lyman a. We therefore explored the influence of the shape of 



1 http://kurucz.harvard.edu/ 



the UV spectrum by studying the impact of the Lyman a line at 
1215 A. We performed a couple of runs with a modified stellar 
spectrum, i.e. a photospheric spectrum with an additional line, 
and the two gas-to-dust ratios (10 and 100). The effects on the 
surface densities of CO, C and C + were found to be negligible. 

3.3. Elemental depletion 

Another possibility to explain the lack of CO and CI in this ob- 
ject is that, even with a normal g/d, the C and O elemental abun- 
dances may be lower than the ones adopted here (e.g. by deple- 
tion onto grains surfaces). We made a model with a stellar UV 
field including UV excess (x = 10 2 ), a standard value of g/d 
(100) and an elemental depletion factor of 10 of heavy elements 
(see Fig. [5] right panel). The amount of CI is in reasonable agree- 
ment with our observation but the CO is largely overestimated. 
The larger amount of CO compared to the g/d=10 and standard 
abundances case is due to mutual shielding by H2. An unlikely 
elemental depletion of ~ 100 is required to match the observa- 
tions of both CO and CI. 

3.4. Spectrum Predictions 

Characteristic flux densities can be recovered by noting that an 
approximate integrated line flux is given by 

S = B v (T ex ,v)n(R olll /D) 2 p6Vcosi 

where B v is the Planck function, T ex is the excitation tempera- 
ture, R out is the outer radius, D is the distance of the star, 5V is 
the local (i.e. sum of thermal + turbulent) line width. For opti- 
cally thin lines p - t, the opa city, and for optically thick lines 
p is of the order of a few (see iGuilloteau & Dutrevl [l998, their 
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q/d rotio=100 UV: y=10 
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Fig. 5. Radial distribution of the surface density of C + , C, and CO for two models which are NOT in agreement with the CO 
observations 



Table 4. Observed and predicted line flux 





CI 3 Pi 




CI S P 2 




Eioo 


P 




(Jy.km/s) 


(Jy.km/s) 


(10 17 


cm -2 ) 


Obs. 


< 6.6 


< 25 




T* 


50 K 


100 K 


50 K 


100 K 






Model A 


30 


48 


68 


135 


12.0 


2.3 


Model B 


13 


15 


33 


52 


2.3 


1.5 


Model C 


26 


37 


60 


111 


8.7 


1.4 


Thick 


60 


128 


125 


300 


1000 


1.5 


Model J / B4 


8- 


11 


24- 


-29 





Model A: x = 10 4 , g/d=100; B: x = 10 2 , g/d=10; C: Kurucz A3, 
g/d=10; Thick: flux for optically thick lines. Predicti ons are for an 
assum ed local line width of 0.2 km/s.Model J is from Jonkheid et al. 
(2007), their model B4. All the line fluxes are scaled for a source dis- 
tance of 140 pc. Eioo and p are the results of the power law fitting of the 
modeled CI surface densities (Fig. [3). 



Eq.2 and Fig.4). In the optically thick case, pdV < 2V out , where 
V out is the Keplerian velocity at the disk outer radius, about 5.5 
km.s -1 in our case. In the optically thin regime, as t oc E/SV, 
where S is the surface density, the emission is just proportional 
to the number of molecules. 

W e used the radiativ e transfer code DISKFIT opt imized for 
disks ( iPietu et al.l 120071: iPavlvuchenkov et alj 120071) . to make 
more accurate predictions for the CI line profiles and integrated 
intensities from our models. The model surface densities were 
extrapolated using simple power laws S(r) = Xioo(//100 AU) _P . 
The disk parameters (size, temperature, Keplerian velocity, incli- 
nation, turbulent width) were taken from Paper I. Our values of 
Eioo and p and the line flux predictions using 7* = 50 and 100 K 
are given in Table|U with the observed upper limits. Tk = 100 K 
is obtained from the best fit model of the J=2-l CO line. 50 K is 
a lower limit, but this temperature provi des a better fit to t he CO 
J=3-2 spectrum obtained at JCMT by iDent et al.l (120051) . Non 
LTE effects are expected to be negligible, because the CI tran- 
sitions have low critical densities of the order o f < 10 4 cm~ 3 
dMonteiro & Flowed [l987t ISchroder et~aD 1 199 lh . while in our 
model, the CI layer is located in a high density region (~ 10 5 
- 10 6 cirT 3 ). 

Considering the uncertainty in the temperature, the non de- 
tection of both CI lines excludes the high CI models (A and C). It 
is just marginally consistent with the low g/d, non photospheric 
UV spectrum model (B), provided the kinetic temperature does 
not exceed 50 K. In this case, the CI lines are mostly optically 
thin, and only the line wings which come from the inner regions 
because of the Keplerian rotation are optically thick. 



3.5. Comparison with other models 

Although there is no specific modeling of the CQ Tau disk and 
star system, it is worth comparing our results with predictions 
from other chemical studies. 

The chemi st ry of HAe disks has also been studied by 
iJonkheid et al.l d2007l) . with different assumptions than in 
our study. The UV trans fer is performed with a 2D code 
(Ivan Zadelhoff et all 120(51 whereas in our case this is a ID 
method. The self-shielding of H2 and CO are calculated assum- 
ing a constant molecular abundance whereas we compute it ex- 
plicitly. Moreover, they mimic the dust growth (and/or settling) 
decreasing the mass of the small grains (i.e., interstellar grains) 
and they consider PAHs. The radiation field is a NEXTGEN 
spectrum (which is a pure photospheric spectrum). Low CO 
amounts can be reproduced by their model B4. They predict that 
CI is the main form of carbon (their Fig. 6), in agreement with 
our model C (photospheric spectrum and low gas-to-dust ratio). 
Jonkheid et al., model B4 yields integrated line intensities of 
about 1-1.3 K.km.s 1 for a beam size of 6.7", an inclination 
of 45° and a turbulent width of 0.2 km.s -1 . This corresponds to 
8 - 1 1 and 24 - 29 Jy.km.s~ 1 for the 492 and 809 GHz lines, 
respectively (Table|4]i. Taking the 1.3 mm flux from CQ Tau and 
using an absorption coefficient of 2 cm 2 .g~' (per gram of dust) 
with a dust temperature of 50 K, the estimated dust mass is 
4.6 10~ 5 M Q . This is likely a lower limit, as the adopted absorp- 
tion coeffici ent is large. The gas mass in the most appropriate 
model from lJonkheid et all (120071) is 10~ 4 M Q , which would im- 
ply a gas-to-dust ratio of ~ 2, only. Our best model (B) with 
g/d=10 predicts somewhat too strong lines. From Table[4] as the 
CI lines are close to the optically thin regime, a reduction of g/d 
by a factor of 2-3 may bring model and observations to agree- 
ment. Thus, whatever chemical model of HAe disk we consider, 
the CO observations and the measured upper limits on both CI 
lines can only be reproduced by a very low gas-to-dust ratio. 

Although a factor of a few may be within the modeling 
uncertainties, it appears difficult to reconcile the current con- 
straints on CO and CI with gas to dust ratio larger than 10. 

The influenc e of X-Rays is not take n into account in our PDR 
code, neither in IJonkheid et al.1 (|2007|) . Although their Lx/Lboi 
are lower, on average, HAe stars a re as strong X-ray emitters as 
TTauri stars dSkinner et al.l 12004). No specific model exist for 
HAe disks, but studies appropriate to T Tauri st ars indicate that 
X Rays are expected to increase the CI intensity (iMeijerink et al.l 
l2008b lGorti & Hollenbachjl2008t lErcolano et al.ll2009l) . We thus 
do not expect the (unknown) X-ray luminosity of CQ Tau to af- 
fect our conclusions. 
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Another alternative would be that models over-predict 
CI due to inaccurate treatment of the physical processes. While 
this cannot be excluded, we note however that the two existing 
HAe disk chemical models make very different approximations 
(one on the radiative transfer, the other on the self-shielding), 
and yet converge towards similar predictions. 

4. Summary 

Our observations of the Cl 3 Pi— > 3 Po and ^Pi^Px transitions 
in the CQ Tau disk leads to significant upper limits which allow 
us to reject chemical models producing a high amount of CI. 
The absence of CI emission together with the low CO emission 
can be explained by currently available chemical models of HAe 
disks, only if the gas-to-dust ratio is notably low (of order 2- 
5) in the CQTau disk. This suggests that the CQTau disk is in 
a transition stage between the proto-planetary and debris disk 
phases. This work also allows us to conclude: 

1 . PDR models which match the CO data are qualitatively com- 
patible with our upper limits, but still overestimate the total 
amount of C I unless the gas-to-dust ratio is of the order of a 
few. 

2. The CI surface density is sensitive to the shape of the UV 
spectrum, with neutral carbon being the dominant species in 
the absence of UV excess. 

3. Introducing the Lyman a line in the stellar UV spectrum does 
not change the CI surface density. 

4. All models matching the low CO and CI intensities predict 
that Carbon should be mainly in the ionized form C + . That 
could be tested by Sofia or Herschel observations. With a 
complete census of CO, CI and C + , the overall surface den- 
sity of gas-phase carbon will be constrained, allowing a bet- 
ter estimate of the g/d ratio. 

While the above conclusions are given under the assumption that 
the chemical models of HAe disks are indeed valid, observing 
more disks, especially large disks with bright CO, would provide 
a useful test of their validity. 
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